NMR relaxation rate in superconducting pnictides: extended s± scenario 
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Recently, several measurements of the nuclear spin lattice relaxation rate T^ ' in the newly discovered super- 
conducting Fe-pnictides have been reported. These measurements generally show no coherence peak below T- 
and indicate a low temperature power law behavior, the characteristics commonly taken as evidence of uncon- 
ventional superconductivity with lines of nodes crossing the Fermi surface. In this work we show that (i) the 
lack of a coherence peak is fully consistent with the previously proposed nodeless extended .s±-wave symmetry 
of the order parameter (whether in the clean or dirty limit) and (ii) the low temperature power law behavior can 
be also explained in the framework of the same model, but requires going beyond the Born model. 
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PACS numbers: 74.20.Rp, 76.60.-k, 74.25.Nf, 71.55.-i 

The recently synthesized^ high-T; superconducting ferrop- 
nictides may be the most enigmatic superconductors discov- 
ered so far One of the biggest mysteries associated with these 
materials is that now, with improved sample quality and sin- 
gle crystal availability, some experiments unambiguously see 
a fully gapped superconducting stat o^'^'^'^'^i'^'^i^i^'' and an s- 
wave pairingii, while others unequivocally point towards line 
nodes in the gapJ^iiiiii^. Particularly disturbing is that on 
both sides data of high quality were reported, by highly rep- 
utable groups, so experimental errors seem unlikely. It is pos- 
sible that full reconciliation will require a highly advanced 
theory that will treat both superconducting and spin density 
wave order parameters on equal footing, and include the in- 
teraction between the two. Nevertheless, it is interesting, and 
important, to investigate more conventional options first. 

Evidence for fully gapped superconductivity comes from 
three different sources: Andreev reflection^i^, exponential 
temperature dependence of the penetration depth a^'^'"^ , and 
the angle-resolved photoemission spectroscopy (ARPES)^*^. 
That three so different probes yield qualitatively the same re- 
sult is very convincing. Yet the nuclear magnetic resonance 
(NMR) spin-lattice relaxation rate, l/Ti, does not show two 
classical fingerprints of a conventional fully-gapped super- 
conductors: the Hebel-Slichter coherence peak and the ex- 
ponential decay at low temperature, but, rather, a power-like 
lav42iii2JdiH^ usually referred to as T^, but in reality some- 
where between and T^^. Such behavior is usually taken to 
be an evidence for a d-wave or similar superconducting state 
with lines of nodes. However, it was pointed ou^ that in 
dirty d-wave samples at low temperatures behavior changes 
from to T (as node lines are washed out into node spots by 
impurities), and this was not observed in ferropnictides. 

So far the evidence in favor of nodeless superconductivity 
seems stronger Therefore, it is interesting to check whether it 
may be possible to explain the results of the NMR experiments 
without involving an order parameter with node lines. 

In this paper we calculate \/T\T for a model superconduc- 
tor consisting of two relatively small semimetallic Fermi sur- 




FIG. 1 : (color online) A depiction of the Fermi surface geometry for 
the Fe-based oxypnictides. Hole and electron FeiTni surface pockets 
are indicated, and an antiferromagnetic wave vector Q is also shown. 



faces, separated by a finite wave vector Q (Fig.[Tl). This is an 
approximation to the Fermi surface of ferropnictides. We in- 
tentionally drop quantitative details that may differ from com- 
pound to compound, and consider the simplest possible case 
with the same densities of states on each surfacaii. We fur- 
ther assume that each surface features the same gapi^, but the 
relative phase between the two order parameters is 7t. This 
is the so-called i+ model, proposed in Ref. [l^ and discussed 
m Refs. [19H20H2II and a number of more recent publications 
by various groups. In the spirit of this model and of model 
calculations we will assume that the total (renormal- 
ized) spin susceptibility is strongly peaked at and around Q. 

We will show here that in this model the Hebel-Slichter 
peak is strongly suppressed already in the clean limit and can 
be entirely eliminated even by a very weak impurity scat- 
tering. On the contrary, the low temperature behavior re- 
mains exponential even in the strong coupling limit. Intro- 
ducing impurities does create strong deviations from the ex- 
ponential behavior, but in the Born approximation the effect 
is stronger just below superconducting temperature and 
weaker at T — > 0, so that the observed power law behavior 
(down to at least 0.1 7"^) is very difficult to reproduce. This 
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behavior, however, can be reproduced if one goes beyond the 
Born Umit of impurity scattering. 

The NMR relaxation rate, assuming a Fermi contact 
hyperfine interaction?-, is given by the standard formula: 
(l/riT) cx limo)^olqlmx±(q,co)/co, where X±{q,&) is the 
analytic continuation of the Fourier transform of the cor- 
relation function x±(r,T) = - ((r,5+(r, -iT)5_ (0,0)) 
averaged (if needed) over the impurity ensemble. Here, 
5±(r, — ix) = exp(//T)5'±(r)exp(— //t), where H is the elec- 
tronic Hamiltonian, T denotes imaginary time, and S± is ex- 
pressed via the electron operators as 5+(r) = \(/j(r)i(/|(r) and 

^-(r) = x|/|(r)x|/|(r). Adopting the above-described model, 
we can keep only the interband contribution to %, in which 
case this formula simplifies to 



lim ImXi2(co)/«, 

0)^0 



(1) 



where Xi2(ki) is obtained by integrating over all q's connect- 
ing the two Fermi surfaces (obviously, only q Q contribute). 
In case of a weakly coupled clean superconductor below Tc we 
have 
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where k and k' lie on the hole and the electron Fermi surfaces, 
respectively, /Sk is the quasiparticle energy in the supercon- 
ducting state, Ai and A2 are the superconducting gaps on hole 
and electron Fermi surfaces, and f{E) is the Fermi distribu- 
tion function. This is a straightforward generalization of the 
textbook expression^^. Following the usual BCS prescription, 
~* EdE/\/E^ — A^, the k-space sum can be converted 
to an energy integral and for a conventional s-wave supercon- 
ductor with Ai = A2 = A one finds 
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The denominator gives rise to a peak just below Tc, the fa- 
mous Hebel-Slichter peak. As pointed out in Ref. [ill it is 
suppressed for the s± state. Indeed, if Ai = — A2 = A, 
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As T decreases from T^, the integral decreases monotonically. 
In a more general case, when Ai A2 < and |Ai | 7^ IA2I, 
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Following Fibich^i^, we assume (— 3/(e)/3e) to be a slow 
varying function and obtain (for Ai > IA2I) 

l/TiT - /(Ai)+2/(Ai,A2)/(Ai)[l-/(Ai)]/r, 
/(Ai,A2) = K(Ai/A2)(Ai+A2)-E(Ai/A2)Ai, 



where K(x) and E(;c) are the complete elliptic integrals of the 
first and second kind, respectively. When Ai = A2, / is re- 
duced to the standard BCS formula, and when Ai = — A2 it 
vanishes identically. 

Let us now include impurity scattering and move to the 
strong coupling limit. Following Samokhin and Mitrovio^^i^, 
we can write down the following formula: 
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For our model gf{(£)) = rii {(o)Zi{(o)(o/Di{(o), gf{(i)) = 

«;(co)(l)/(co)/D,-(co), where A-(co) = ^[Zi(co)co]^-(j)?((o), 
Z,((o) is the mass renormalization, (|),(co) = Z,(a))A,(a)) and 
«/ (co) is a partial density of states. 

The renormalization function Z, (co) and complex order pa- 
rameter (|),(co) have to be obtained by a numerical solution 
of the Eliashberg equations. On the real frequency axis they 
have the form (we neglect all instant contributions and con- 
sider a uniform impurity scattering with the impurity potential 
Vij = v). 
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where ® = 1 -a+a{ [gf(co) +gf («)]' - [gf{(o)~g^i(o)f}, 

Y= 2co/kN{0) is the normal state scattering rate, A^(0) is 
the the density of states at the Fermi level, c is the impurity 
concentration, and o = [7cw(0)t] impurity strength 

l + [jriV(0)v]- f J fe 

(0^0 corresponds to the Born limit, while o = 1 to the 
unitary one). Kernels Kjj{z,(i)) are 
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where Bij{Q.) is equal to Bij{Q.) in the equation for A, and to 
\Bjj{Q.) I in the equation for Z. Note that all retarded interac- 
tions enter the equations for the renormalization factor Z with 
a positive sign. 

It is well known that pair-breaking impurity scattering dra- 
matically increases the subgap density of states just below Tc, 
and even weak magnetic scattering can eliminate the Hebel- 
SUchter peak in conventional superconductors. In our model, 
the same effect is present due to the nonmagnetic interband 
scattering (the magnetic scattering, on the contrary, is not pair- 
breaking in the Born limit)2^. Since the Hebel-Slichter peak 
is not present in this scenario even in a clean sample, the pair- 
breaking effect is more subtle: it changes exponential behav- 
ior below Tc to a more power-law like one (the actual power 
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and extent of the temperature range with a power-law behav- 
ior depend on the scattering strength). Note that in the Born 
approximation the exponential behavior is always restored at 
low enough temperature, unless the impurity concentration is 
so strong that Tc is suppressed by at least a factor of two^. 

Another well-known pair-breaking effect is scattering by 
thermally excited phonons (or other bosons). This is, of 
course, a strong coupling effect. For instance, strong cou- 
pling can nearly entirely eliminate the Hebel-Slichter peak in 
a conventional superconducto r^^i^*^ . However, this effect is 
even more attached to a temperature range just below Tc since 
at low temperature boson excitations are exponentially sup- 
pressed. 

Currently, most experimental data for ferropnictides go 
down in temperature to ^ 0.2 — Q3Tc, but some results are 
available at temperatures as low as 0.1 r^. So far exponential 
behavior has not been observed, which casts doubt that Born 
impurity scattering may be responsible for such behavior 

Unitary scattering, on the other hand, has rather different 
low-temperature behavior As discussed, for instance, by Pre- 
osti and Muzikar-", unitary scattering in the case of s± su- 
perconducting gap (our choice of |Ai| = IA2I corresponds to 
their parameter r set to zero) the subgap density of states is 
controlled by the unitarity parameter o, while the suppression 
of Tc is controlled by a different parameter, namely, by the 
net scattering rate y. The unitary limit corresponds to o ^ 1, 
but Y may be rather small at low concentrations. The physi- 
cal meaning is that here the diluted unitary limit corresponds 
to the so-called "swiss cheese" model: each impurity creates a 
bound state that contributes to the subgap density of states, but 
hardly to the Tc suppression. Indeed, Preosti and Mizikar have 
shown (see Fig. 1 in Ref.lJl}) that in this limit nonzero density 
of states at the Fermi level appears already at zero tempera- 
ture at arbitrary low impurity concentration. That is to say, 
the bound state has zero energy. This is a qualitatively new 
effect compared to the Born limit: In a dilute unitary regime 
(Y <C A), the NMR relaxation at T « is mainly due to the 
bound states at £ = 0; upon heating \/T\T initially remains 
constant, or may even slightly decrease because if depopulat- 
ing of the bound state. When the temperature increases fur- 
ther, at some crossover temperature the relaxation becomes 
dominated by thermal excitations across the gap and \/T\T 
starts growing exponentially. When the gap is suppressed by 
the temperature as to become comparable with Y, yet another 
effect kicks in, broadening of the coherence peak near Tc (less 
important for our s± state). Thus, unitary scattering makes 
the \/T\T temperature dependence rather complex, although 
the strongly unitary regime with low impurity concentration 
is rather far from a power law behavior (even though being 
strongly nonexponential, see Fig.|3]l. 

These qualitative arguments suggest that neither purely 
Born nor purely unitary limits are well suited for explain- 
ing the observed \/T\ behavior: the former leads to an ex- 
ponential behavior at low temperatures, the latter to Kor- 
ringa behavior On the other hand, an intermediate regime 
seems to be rather promising in this aspect. Indeed, the en- 
ergy of the above-mentioned bound state is related to a as 
Eh = |A|\/1 a. Thus, by shifting o toward an intermedi- 




FIG. 2: (color online) Temperature dependence of the spin-lattice 
relaxation rate l/T^T calculated in the clean limit (y = 0) for a con- 
ventional .s-wave superconductor (Ai = A2) and for an s± supercon- 
ductor (Ai = — A2). Results of both weak and strong coupling ap- 
proximations are shown. 



ate scattering o ^ 0.5 (which is probably more realistic than 
either limit anyway) and increasing y we create broad and 
monotonously increasing density of states, much closer to 
what would be expected from the NMR data. It is worth not- 
ing that a distribution of a's (presence of different impurities 
with different scattering strength), which we do not consider 
here, will also lead to broadening of the bound state and work 
as an effectively enhanced y (therefore it is reasonable to try 
relatively large y's keeping in mind that part of this is simulat- 
ing). 

We will now illustrate the above discussion using spe- 
cific numerical models. First, we present numerical solu- 
tions of the Eliashberg equations using a spin-fluctuation 
model for the spectral function of the intermediate boson: 
B,/(co) = 'kijTzQ.sf/i'^lf + CO^), with the parameters Q.sf = 25 
meV, A,ii = A,22 = 0.5, and X12 = A,2i = —2. This set gives a 
reasonable value for Tc ~ 26.7 K. A similar model was used 
in Refill to describe optical properties of ferropnictides. The 
actual details of the function are in fact not important; our us- 
age of this particular function does not constitute an endorse- 
ment or preference compared to other possibilities, but is just 
used here for concreteness. 

In Fig.|2]we compare the temperature dependence of the re- 
laxation rate calculated as described above in the clean limit 
for a conventional s-wave superconductor (Ai = A2) and for 
an s± superconductor (Ai = — A2), both in the weak and in 
the strong coupling Umits. We observe that while in the con- 
ventional case strong coupling makes a big difference by sup- 
pressing the coherence peak, in the s± state, where no coher- 
ence effects take place, strong coupling is not really impor- 
tant. In Fig.[3ja) we show the effect of impurities in the Born 
limit. We have found that for an impurity scattering of the 
order of OATcq, where Tcq is the transition temperature in the 
absence of impurities, there is a moderate suppression of Tc 
(less than 20%). More importantly, the strong deviation from 
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FIG. 3: (color online) Calculated temperature dependence of the 
spin-lattice relaxation rate \/T\T for an s± superconductor in the 
strong coupling approximation without impurities (y = 0), with im- 
purities in the Born limit (y = CSr^), and in the unitary limit for 
small (y= 0.07Ao) and large (y = 0.35Ao) impurity concentrations. 
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FIG. 4: (color online) Calculated spin-lattice relaxation rate for the 
s± superconducting state together with experimental l/Ti data from 
several groups, as indicated. T--^ trend is also shown. 



exponential behavior in 1 /T\T appeared. Above ~ 0.27"^ the 
dependence can be well represented by a power law, but with 
an exponent closer to 5.5 for clean, and 4 for dirty samples 
(as opposed to the experimentally observed 1.5 to 2). Further 
increase of the scattering rate leads to a too strong suppression 
of transition temperature and a too large relaxation rate right 
below Tc- Thus, impurity scattering in the Born limit can not 
fully explain the NMR data. 

As shown in Fig.[3lb), the unitary limit also does not repro- 
duce the experimental data (represented by the approximate 
behavior \/T\T °^ T^) neither for small nor for large impurity 



concentrations, and even predicts a slight non-monotonicity 
for small y = O.OTAq. Here, Aq is the low temperature value 
of the energy gap without impurity scattering. 

On the other hand, the experimental results can be re- 
produced very well if one assumes the intermediate regime 
of impurity scattering. Fig. 2] shows various experimental 
dat a'^''^''^ together with our calculations for the s± gap with 
a taken as 0.4 and interband y taken as O.SAq. This y corre- 
sponds to a relatively dirty superconductivity, but the effect 
of interband scattering on Tc, for given y, in this s± state is 
smaller than would be effected by intraband scattering of the 
same magnitud o^^i^^ . Besides, as mentioned above, a distri- 
bution of o's will lead to a similar broadening of the DOS 
for smaller y. We observe again that the s± state exhibits no 
coherence peak. As opposed to the Born and unitary limits, 
intermediate-a scattering is capable of reproducing the exper- 
imental behavior, usually described as cubic, but in fact prob- 
ably closer to T^-^ in 1 /T\. Note that there is no universality 
after the 2.5 power of T, it is simply result of a fitting. 

We want to emphasize that this analysis does not prove 
that the origin of the power law behavior is dirty-limit 
intermediate-a scattering in an s± state. It is fairly possible 
that more complex physics, possibly related to coexistence of 
superconductivity and spin density wave order, plays a role. 
But it clearly demonstrates that such a behavior does not prove 
existence of gap nodes on the Fermi surface. 

To summarize, we have shown that the lack of a coherence 
peak is very naturally explained in the framework of the s± 
superconducting state even in the clean limit, and even more 
so in the presence of impurities. However, a clean s± su- 
perconductor would show an exponential decay of the relax- 
ation rate 1 /T\ below 7^, contrary to what has been observed 
in NMR experiments. Strong coupling effects and impurity 
scattering in Born approximation transform this exponential 
behavior into a power law-like for temperatures T > Q.2Tc, 
but it is difficult to reproduce the actual experimental temper- 
ature dependence. On the other hand, an intermediate-limit 
scattering (neither Born nor unitary), can reproduce experi- 
mental observations rather closely. While we did not address 
in this paper any effects that strong scattering may have on the 
other physical properties (this is left for future publications), 
we want to emphasize that there is an important difference be- 
tween the scattering effect on the properties related to q = 
response (penetration depth, tunneling, specific heat) and l/Ti 
that probes mainly the q ~ Q response. 

We would like to thank A. Chubukov, I. Eremin, K. Ishida, 
and G.-q. Zheng for useful discussions. 

Note added As this paper was being finalized, we became 
aware of a work by Chubukov et alM- who arrived at similar 
conclusions using a different approach. After submission of 
our paper. Bang and Choi^^ reported a similar, but indepen- 
dent research, again reaching the same conclusion as ours. 
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